Mediation of Long Range Charge Transfer by Kondo Bound States 
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We present a theory of non-equilibrium long range charge transfer between donor and acceptor cen- 
ters in a model polymer mediated by magnetic exciton (Kondo) bound states. Our model produces 
electron tunneling lengths easily exceeding lOA, as observed recently in DNA and organic charge 
transfer systems. This long ranged tunneling is effective for weak to intermediate donor-bridge 
coupling, and is enhanced both by weak to intermediate strength Coulomb hole-electron attraction 
(through the orthogonality catastrophe) and by coupling to local vibrational modes. 



PACS Indices: 72.15.Qm, 85.65.-|-h,87.15.-v 

Charge transfer (CT) via quantum mechanical tunnel- 
ing of electrons or holes through a molecular bridge be- 
tween well separated donor (D) and acceptor (A) sites, 
is central to many biological processes [l]-H, and may 
play a key role in molecular electronics applications p,p| . 
For proteins, the CT rate typically decays exponentially 
in the D/A separation with the inverse tunneling length 
f3 = 1.1 ± 0.2A~^ H. However, evidence for anoma- 
lously small (3 values ( of order O.lA^^) has accrued for 
CT between metallic intercalating complexes or organic 
complexes in DNA [pM, CT from ferrocene attached to 
conjugated polymers ^^, and a few other systems 0. 
Values of /3 ~ 1A~^ are well accounted for in effective 
one-electron theories with the D/A levels typically resid- 
ing in the middle of the large insulating (4-6 eV) gaps of 
the polypeptide bridges of proteins 1^. Understanding 
smaller (3 values has generally required "fine tuning" of 
the D/A levels close to the lowest unoccupied molecu- 
lar level (LUMO) [§. For DNA, this issue is coupled to 
the magnitude of the electronic energy gap in DNA, an 
exciting and hotly debated subject [[l0| - [l3| . 

In this paper, we show that long range CT, with tun- 
neling lenghths (3~^ > 10 A can be achieved in model 
systems through mediation by a semiconducting, molec- 
ular analogue of the Kondo state in which a localized 
electron or hole moment is antiferromagnetically cou- 
pled to a hole or electron donated to the bridge p4[ . 
While our model for the tunneling bridge is a conju- 
gated polymer, the conclusions are generic for a gapped 
molecule. The bound state mechanism provides a ro- 
bust basis for "fine tuning" D levels close to the LUMO. 
We carry out our calculations both with a semi-analytic 
variational method jlq] as well as the density matrix 
renormalization group (DMRG) approach, and the re- 
sults agree very well. Adding the dielectrically screened 
electron-hole Coulomb interaction to the Kondo coupling 
leads to an enhancement of tunneling for weak to in- 
termediate Coulomb strength through the orthogonality 
catastrophe, and a crossover to deep Coulomb binding 
for larger Coulomb interaction. Coupling to local vibra- 
tional modes of the metal complexes also enhances the 
CT rate, and can lead to normal or inverted regimes de- 



pending upon the strength of the binding and the D/A 
distance. 

CT via a single tunneling event is intiated out of equi- 
librium either by the absorption of light (photolysis) or 
transfer of charge by radicals (radiolysis) so that an elec- 
tron or hole at the D complex is energetically excited 
above an empty level in the A complex; in this paper we 
focus on photolytically prepared charge transfer. In the 
simplest "diabatic" theory for which the electron tun- 
neling matrix element Hoa is small, the CT rate k^ is 
well described by a golden rule/linear response approach 
(Marcus theory |l|), with 

kctiT)^^\HDA\'FiT) (1) 

where F{T) is the Franck-Condon factor M. Hjja, which 
is the main subject of this paper, carries the dominant 
distance dependence. 

The CT scenario involves three steps, and can only 
arise when, as in Fig. 1, the excited D complex energy 
lies at or near the bridge LUMO (a similar situation can 
be realized for hole transfer with suitably reversed ener- 
getics). Upon excitation (step 1), the D complex donates 
an electron to the bridge, leaving behind a remnant hole 
spin with antiferromagnetic coupling to the bridge elec- 
tron (this interaction drives the Kondo binding in step 
2). In step 3, the electron is transferred to the distant 
donor level. 

Experimental realization of the appropriate D energet- 
ics seems possible. For example, for the Ru-complex met- 
allointercalator used in DNA CT experiments g], a 480 
nm photoexcitation transfers electrons from Ru(II) to 
the ligand (L*) from which electron transfer to the DNA 
presumably occurs (there is no direct Ru-DNA overlap), 
leaving a hole spin in the t2g Ru shell. For DNA gaps of 
1-2 eV, it is plausible that the excited Ru(III)-L* com- 
plex lies above the electron addition threshold, though 
estimates of the gap range from eV |1^JI^, 1-2 eV 
0,0, and on up to 7 eV [0. The 3.9 eV photoexcita- 
tion energy of the adenine isomer 2-aminopurine (Ap) is 
very comparable or in excess of the HOMO-LUMO gap 
for isolated bases [0, although the initial CT matrix el- 
ement to neighboring DNA bases is likely weaker than 



in the metallointercalator case, for which the ligand re- 
sides between adjacent base pairs with more effective tt- 
overlap. Excitation of Rh(phi)2(bpy)'^^ will produce con- 
ditions favorable for hole transfer to adjacent sequences 
of DNA-A bases @|9l. For the 5-(p-tolyl)-tetracene 
donor of ref. B , the D photoexcitation energy is nearly 
equal to the p-phenylenevinylene bridge gap for bridges 
of 3-5 units, and judging from the trend, should become 
higher for longer bridges. Thus, our considerations are 
clearly relevant to these systems. 

With these energetics, the Kondo singlet bound state 
(KS) [p^ is the most stable charge zero singlet state, 
though the charge ±1 spin 1/2 state (free moment) can 
be globally stable in a molecule or semiconductor pl| . 
Rapid KS formation is required to compete with other 
processes, as discussed in Fig. 1. A lower bound estimate 
for KS formation is tk ~ h/Es , where Eb is the binding 
energy, in analogy to the reasoning used for Kondo reso- 
nance formation time in quantum dots p2| . tk is of the 
order of a few femtoseconds for Eb = 0.1 — 0.2eV. In the 
semiconducting/molecular context, it will be divided by 
a Franck-Condon factor Fk since the excited D level re- 
sides above the bound state. Provided Fk < 0.001 p4[ , 
KS formation will likely beat luminescence (nanosecond 
time scale) and phonon radiative decay to the A state 
(large polaron formation) which can be estimated, utiliz- 
ing polyacetylene parameters [g3|, to be of order 100 ps. 
We stress that (i) no simple polaron dynamics will cause 
kct to decay exponentially in the D/A separation, and 
(ii) any polaron mediated contribution to kd is bounded 
above by the large polaron result |2^ . 

We simply model the insulating bridge as a dimerized 
tight binding chain of length N {N even) with Hamilto- 
nian He given by p3] 
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N 

Y,{t+i-iyt'){cl,^c,, + h.c.) (2) 



where cl_^_i^ creates an electron of spin a at site i. This 
model describes the tt^ bands of polyacetylene, and can 
easily be extended to other conjugated polymers (as in 
refs. ||J§1), but the bound state formation is not special 
to this model bridge. For large N, this model produces 
two bands with a full width of 4t and a gap of 2A = 4t' . 
The acceptor complex is placed adjacent to site N of 
the chain, and modeled as an Anderson impurity in the 
empty orbital regime [l3, with energy ea- To simplify 
the analysis, we have taken the Coulomb repulsion U 
between D holes to be infinite, and as in the Ru case 
we assume there is a single hole in an otherwise filled 
shell, with energy relative to the filled shell of e/j. The 
hybridization Hamiltonian coupling the D/A complexes 
and bridge is 



where \DQ > is the filled shell state of the Donor and 
\Da > the state of a single hole and Xq^- = \D0 >< Dcr\. 
We choose en above the LUMO, and ca to be within the 
gap. 

We compute Hjja by first identifying the appropriate 
two state system from diagonalizing with either: (1) the 
DMRG approach, which is exact to within targeted nu- 
merical precision pfl], and (2) a version of the variational 
wave function approach for magnetic impurities in met- 
als pioneered by Varma and Yafet as well as Gunnarsson 
and Schonhammcr |13]. We build our variational states 
by first diagonalizing the bridge with the non-interacting 
A "impurity" included, and then employ the Ansatz: 



JO- 



a>]\BA> 



(4) 



Ti-BDA = ^[Vo-'^OcrCl.-CT + VAalcNa + h.c] 



(3) 



where bj^ creates an electron in orbital j of the bridge-A 
system, aj are variational coefficients, and the particle- 
hole excitations are controlled both by the spin degen- 
eracy Ns (here Ng — 2) and the bridge gap. As is 
well known Mm , each added particle-hole excitation con- 
tributes a term of relative order l/N^ to the energy. This 
variational ansatz does an excellent job of describing the 
Kondo state in metals (the second term captures the 
screening of the local moment, and as such is nonper- 
turbative as opposed to, e.g., direct expansion in U or 
MoUer-Plesset methods). Because of the Hilbert-Space 
truncation associated with large U, the above wave func- 
tion cannot be reduced to a single determinant by a uni- 
tary transformation. Hence, for our model, DMRG gives 
an "exact configuration interaction" solution while the 
variational approach provides a physically motivated ba- 
sis set reduction solution. For the purposes of calculating 
HdA: the methods agree exceedingly well. 

Given the two lowest "adiabiatic" eigenstates, we next 
identify the out-of-equilibrium "diabatic" states corre- 
sponding to the excited electron predominantly localized 
to the left (D) side, or right (A) side. The Generalized 
Mulliken-Hush method [2^ accomplishes this by finding, 
within the two-state space, the eigenstates for the elec- 
tric dipole operator, which maximally localize the charge 
to left or right. Hoa is then approximated by the off- 
diagonal matrix element of the two-state Hamiltonian in 
the dipole-diagonal basis. 

Our results for H^a vs. bridge length N are shown in 
Fig. 2(a), where we have set Vd = Va for convenience, 
and chosen eu well above the LUMO, and Vd is of order t 
p8[ . We note the following features of our curves: (i) for 
short distances, Hda scales with VdVa, perturbatively, 
so that stronger initial coupling yields stronger tunneling. 
On the other hand, for large distances, the weaker Vd 
gives better tunneling, because the more loosely bound 
KS cloud has better overlap on the distant A complex, 
(ii) By taking the bridge spacing to be of the order of SA 
as in DNA, the slope for our curve with Vd = 0.75t at 



long distances corresponds to (3 ~ O.lA^^, and this corre- 
sponds to a binding energy of the order of 0.05-0.2 eV (see 
Fig. 3) . Obviously our model is too simple to realistically 
describe CT in DNA; rather, we stess the ready achiev- 
ment of small /3 values, (iii) The agreement between the 
DMRG calculations (x's) and variational (lines) is quanti- 
tatively excellent in the long distance regime where Hjjji 
is quite small, and good in the short distance regime. 
Hence, our variational approach is a very useful, accu- 
rate, physically motivated basis set reduction scheme for 
this model, (iv) By varying i', Vd, cd we can easily ob- 
tain the typical protein results (/3 ~ 1A~^) in our model. 
Analytic solutions are possible in two limits. First, in a 
continuum limit in which the D/A complexes are placed 
at separation Na in a one-dimensional wire, and Eb << 



E 



^i^ exp{-^2m*EBNa/h), 
H'l\2{t^ - i'2)a2] and the 



t', we find that Hjja 
with effective mass m* = fi^t' /[2{t^ 
Kondo binding energy Eb = Vj^/(eo|e£) — 2t'p), with 
eo — h /{%m*a^). The prefactor reflects a "slave bo- 
son" renormalization of the matrix element Vd pS] , and 
the exponent arises from A/KS overlap. Thus there is a 
maximum Hba for a given bridge length, though this is 
practically relevant for only long bridges given the weak 
prefactor dependence upon Eb- Second, for a discrete 
bridge Hua is expressed approximately as 



\H 



DA\ 



WdVaY. 



ek-2t' + EB ' 



(5) 



where ip{k,i) is the wave function amplitude at site i 
of a bridge-only conduction orbital of energy efc (no A 
complex). This describes the propagation from the D 
complex to bridge site 1, then to bridge site N and onto 
the A complex (hole transfer processes are also possible 
but less dominant). 

Given the semiconducting character of the bridge, the 
dielectrically screened Coulomb interaction between the 
D hole and the bridge electron is of manifest importance. 
This is modeled by the interaction 



He 



N+l ^2 






(6) 



where X„„ = \Da >< D<j\, and n^ is the electron oc- 
cupancy at site i, and e is the high frequency dielectric 
constant. There are two important effects: (i) For suffi- 
ciently small interaction {Vc = e^/(ea) < 1 eV in prac- 
tice) , the modification of the bridge wave functions leads 
to a multiplication of Vd by < BA\BA > where \BA > is 
the filled valence sea in the presence of the D hole. This 
is directly analogous to the orthogonality catastrophe of 
the core level x-ray absorption problem. The effects of 
this on Hda are shown in Fig. 2(b). Effectively, small 
Vc reduces Vd and thus enhances Hda- (ii) For larger 
Vc, Coulomb binding overtakes Kondo coupling, and the 
stronger binding reduces Hda- 



Finally, we introduce coupling to local vibrational 
modes of the D/A complexes. We make the harmonic 
approximation for the phonons, and treat the displace- 
ments classically, adding to our model 



Hel,ph = X^^ixDXaa + XAaj^tta) 



(7) 



where A is the electron-phonon coupling constant, xd 
and xa are the nuclear coordinates, and we add a har- 



monic potential 4- (x 



^l 



After the substitu- 



tion xd/a -^ X/Kxd/Ai with dimensionless coordinates 
xd/Ai the Hamiltonian depends only on a single parame- 
ter }? jK- By neglecting electron-phonon coupling within 
the bridge [gO| , we have calculated the total electronic en- 
ergy of the groundstate and first excited state (two state 
model) as a function of the nuclear coordinates. These 
constitute the Born-Oppenheimer potential energy sur- 
faces (PES) for the nuclear motion of the two Ander- 
son impurities. The resonance splitting at the transition 
state is 2Hda %- 

In Fig. 3, we compare with the electron-only MuUiken- 
Hush approach. Like the Coulomb interaction discussed 
previously, coupling to phonons leads to an enhanced CT 
rate. The MuUiken-Hush result corresponds to the limit 
of vanishing phonon coupling. In the perturbative regime 
y? jK <K 1, we find a hnear dependence of ln{HDA/t) on 
X^/K. We can roughly explain this by replacing ec in 
the asymptotic, electron-only problem with renormalized 
energy level eu = cd + X^ /Kxd and expanding In Hda 
to linear order in X^ /K. This result varies little with xd, 
since the reaction path approximately follows xa- 
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FIG. 1. Dynamics of Kondo bound state (KS) mediated 
CT. Photolysis leads to the excited state D" from which an 
electron is transferred to the bridge leaving a remnant D 
hole spin. For sufficiently rapid KS formation, KS mediated 
CT (1) is more rapid and likely than luminescence (2) or 
phonon-radiative decay to the acceptor (3). 
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FIG. 2. Bridge length (A'^) dependence of electron tun- 
nelling matrix element Hda (a) Dependence of Hda on 
bridge-donor coupfing. Lines; variational method. Symbols: 
DMRG approach. Dotted line: to and Vd have been reduced 
so as to yield the same continuum limit Kondo binding en- 
ergy as the solid line. (B) Effect of Coulomb interaction Vc 
between the bridge electron and remnant D hole. 
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FIG. 3. Distance dependence of Hda'- results from the 
splitting of the Born-Oppenheimer surfaces are compared 
with the Mulliken-Hush algorithm as a function of base pairs 
N for parameters eo ~ 2f, ea ~ 0.5t, V — 0.75t. Thick solid 
lines: exponential decay assuming a — SAA 



